Characterisation and representation of non-dissipative electromagnetic medium with a double light cone

We study Maxwell's equations on a 4-manifold N with a medium that is non-dissipative and has a linear and pointwise response. In this setting, the medium can be represented by a suitable (2,2)-tensor on the 4-manifold N. Moreover, in each cotangent space on N, the medium defines a Fresnel surface. Essentially, the Fresnel surface is a tensorial analogue of the dispersion equation that describes the response of the medium for signals in the geometric optics limit. For example, in isotropic medium the Fresnel surface is at each point a Lorentz light cone. In a recent paper, I. Lindell, A. Favaro and L. Bergamin introduced a condition that constrains the polarisation for plane waves. In this paper we show (under suitable assumptions) that a slight strengthening of this condition gives a pointwise characterisation of all medium tensors for which the Fresnel surface is the union of two distinct Lorentz null cones. This is for example the behaviour of uniaxial medium like calcite. Moreover, using the representation formulas from Lindell et al. we obtain a closed form representation formula that pointwise parameterises all medium tensors for which the Fresnel surface is the union of two distinct Lorentz null cones. Both the characterisation and the representation formula are tensorial and do not depend on local coordinates

Drifting instabilities of cavity solitons in vertical cavity surface-emitting lasers with frequency selective feedback

In this paper we study the formation and dynamics of self-propelled cavity solitons (CSs) in a model for vertical cavity surface-emitting lasers (VCSELs) subjected to external frequency selective feedback (FSF), and build their bifurcation diagram for the case where carrier dynamics is eliminated. For low pump currents, we find that they emerge from the modulational instability point of the trivial solution, where traveling waves with a critical wavenumber are formed. For large currents, the branch of self-propelled solitons merges with the branch of resting solitons via a pitchfork bifurcation. We also show that a feedback phase variation of 2\pi can transform a CS (whether resting or moving) into a different one associated to an adjacent longitudinal external cavity mode. Finally, we investigate the influence of the carrier dynamics, relevant for VCSELs. We find and analyze qualitative changes in the stability properties of resting CSs when increasing the carrier relaxation time. In addition to a drifting instability of resting CSs, a new kind of instability appears for certain ranges of carrier lifetime, leading to a swinging motion of the CS center position. Furthermore, for carrier relaxation times typical of VCSELs the system can display multistability of CSs.Comment: 11 pages, 12 figure

Flavor composition of neutrinos from choked gamma-ray bursts

Choked gamma-ray bursts (CGRBs) are possible neutrino sources that have been proposed as capable of generating the flux detected by IceCube, since no accompanying gamma-ray signal is expected, as required by observations. We focus on obtaining the neutrino flux and flavor composition corresponding to CGRBs under different assumptions for the target photon density and the magnetic field of the emission region. We consider the injection of both electrons and protons into the internal shocks of CGRBs, and using a steady-state transport equation, we account for all the relevant cooling processes. In particular, we include the usually adopted background of soft photons, which is the thermalized emission originated at the shocked jet head. Additionally, we consider the synchrotron photons emitted by the electrons co-accelerated with the protons at the internal shocks in the jet. We also obtain the distributions of pions, kaons, and muons using the transport equation to account for the cooling effects due not only to synchrotron emission but also interactions with the soft photons in the ambient. We integrate the total diffuse flux of neutrinos of different flavors and compute the flavor ratios to be observed on Earth. As a consequence of the losses suffered mainly by pions and muons, we find these ratios to depend on the energy: for energies above ~(10^5-10^6) GeV (depending on the magnetic field, proton-to-electron ratio, and jet power), we find that the electron flavor ratio decreases and the muon flavor ratio increases, while the tau flavor ratio increases only moderately. Our results are sensitive to the mentioned key physical parameters of the emitting region of CGRBs. Hence, the obtained flavor ratios are to be contrasted with cumulative data from ongoing and future neutrino instruments in order to assess the contribution of these sources to the diffuse flux of astrophysical neutrinos.Comment: 23 pages, 18 figures. Accepted for publication in A&

New Perspective on Galaxy Clustering as a Cosmological Probe: General Relativistic Effects

We present a general relativistic description of galaxy clustering in a FLRW universe. The observed redshift and position of galaxies are affected by the matter fluctuations and the gravity waves between the source galaxies and the observer, and the volume element constructed by using the observables differs from the physical volume occupied by the observed galaxies. Therefore, the observed galaxy fluctuation field contains additional contributions arising from the distortion in observable quantities and these include tensor contributions as well as numerous scalar contributions. We generalize the linear bias approximation to relate the observed galaxy fluctuation field to the underlying matter distribution in a gauge-invariant way. Our full formalism is essential for the consistency of theoretical predictions. As our first application, we compute the angular auto correlation of large-scale structure and its cross correlation with CMB temperature anisotropies. We comment on the possibility of detecting primordial gravity waves using galaxy clustering and discuss further applications of our formalism.Comment: 10 pages, 2 figures, accepted for publication in Physical Review

Current cosmological constraints from a 10 parameter CMB analysis

We compute the constraints on a ``standard'' 10 parameter cold dark matter (CDM) model from the most recent CMB and data and other observations, exploring 30 million discrete models and two continuous parameters. Our parameters are the densities of CDM, baryons, neutrinos, vacuum energy and curvature, the reionization optical depth, and the normalization and tilt for both scalar and tensor fluctuations. Our strongest constraints are on spatial curvature, -0.24 < Omega_k < 0.38, and CDM density, h^2 Omega_cdm <0.3, both at 95%. Including SN 1a constraints gives a positive cosmological constant at high significance. We explore the robustness of our results to various assumptions. We find that three different data subsets give qualitatively consistent constraints. Some of the technical issues that have the largest impact are the inclusion of calibration errors, closed models, gravity waves, reionization, nucleosynthesis constraints and 10-dimensional likelihood interpolation.Comment: Replaced to match published ApJ version. More details added. 13 ApJ pages. CMB movies and color figs at http://www.hep.upenn.edu/~max/10par_frames.html or from [email protected]

Accelerating solid-fluid interaction using Lattice-Boltzmann and Immersed Boundary coupled simulations on heterogeneous platforms

We propose a numerical approach based on the Lattice-Boltzmann (LBM) and Immersed Boundary (IB) methods to tackle the problem of the interaction of solids with an incompressible fluid flow. The proposed method uses a Cartesian uniform grid that incorporates both the fluid and the solid domain. This is a very optimum and novel method to solve this problem and is a growing research topic in Computational Fluid Dynamics. We explain in detail the parallelization of the whole method on both GPUs and an heterogeneous GPU-Multicore platform and describe different optimizations, focusing on memory management and CPU-GPU communication. Our performance evaluation consists of a series of numerical experiments that simulate situations of industrial and research interest. Based on these tests, we have shown that the baseline LBM implementation achieves satisfactory results on GPUs. Unfortunately, when coupling LBM and IB methods on GPUs, the overheads of IB degrade the overall performance. As an alternative we have explored an heterogeneous implementation that is able to hide such overheads and allows us to exploit both Multicore and GPU resources in a cooperative way